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The collective excitation phenomena in atomic nuclei are studied in two different 
formulations of the Random Phase Approximation (RPA): (i) RPA based on cor- 
related realistic nucleon-nucleon interactions constructed within the Unitary Corre- 
lation Operator Method (UCOM), and (ii) relativistic RPA (RRPA) derived from 
effective Lagrangians with density-dependent meson-exchange interactions. The for- 
mer includes the dominant interaction-induced short-range central and tensor cor- 
relations by means of an unitary transformation. It is shown that UCOM- RPA 
correlations induced by collective nuclear vibrations recover a part of the residual 
long-range correlations that are not explicitly included in the UCOM Hartree-Fock 
ground state. Both RPA models are employed in studies of the isoscalar monopole 
resonance (ISGMR) in closed-shell nuclei across the nuclide chart, with an emphasis 
on the sensitivity of its properties on the constraints for the range of the UCOM 
correlation functions. Within the Relativistic Quasiparticle RPA (RQRPA) based 
on Relativistic Hartree-Bogoliubov model, the occurrence of pronounced low-lying 
dipole excitations is predicted in nuclei towards the proton drip-line. From the anal- 
ysis of the transition densities and the structure of the RQRPA amplitudes, it is 
shown that these states correspond to the proton pygmy dipole resonance. 
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I. INTRODUCTION 

Among various theoretical approaches to nuclear structure, two pathways have been 
extensively pursued over the past decades: (i) models based on effective nuclear interactions 
constrained by the properties of nuclear matter and bulk properties of finite nuclei (e.g. 
Skyrme [lj], Gogny [2], and relativistic models based on exchange of effective mesons Q), 
and (ii) models which start from a realistic nucleon-nucleon (NN) interaction. Recently, 
several modern realistic NN interactions have been constructed, e.g., the Argonne V18 15], 

n n 

the CD-Bonn |4|], and chiral potentials p, which reproduce the experimental NN phase- 
shifts with high accuracy. Within ab initio Green's function Monte Carlo 17] and no-core 

n 

shell model calculations [a] of ground state properties and low-lying excitation spectra of 
light nuclei it was shown that realistic NN interactions, supplemented by a three-nucleon 
force, allow for a quantitative description of experimental data l| . 

Realistic NN interactions cannot be directly employed in a standard Hartree-Fock (HF) 
scheme due to the importance of interaction-induced correlations in the many-body state 
beyond the simple HF Slater determinant. Therefore, an effective, phase-shift equivalent 
interaction has to be derived by explicitly accounting for the dominant correlations. One way 



to tackle this issue is the Unitary Correlation Operator Method (UCOM) which descri 



Des the 



11 



12, 



13. 



short-range central and tensor correlations by means of an unitary transformation 
Q]. The unitary transformation of the Hamiltonian including a realistic NN potential results 
in a correlated effective interaction well suited for the application with simple uncorrelated 
many-body states. An alternative method to derive phase-shift equivalent, low- momentum 
effective interactions is the V[ ow -fc renormalization group approach [l5^ . 

Studies of collective excitation phenomena in atomic nuclei provide valuable insight into 
many properties of the underlying effective interactions employed in solving the nuclear 
many-body problem. In order to describe small-amplitude collective excitations within the 
UCOM framework, one can employ the random-phase approximation (RPA) Q], based on 
the HF single-nucleon basis. Since three-body interactions presently are not included, the 
results provide information on their importance for the understanding of collective nuclear 
excitations. The UCOM-RPA model can also be employed to evaluate the contributions of 
RPA correlations to the ground state energy which go beyond the mean-field picture H. In 
addition to the short-range correlations described explicitly by the unitary transformation, 
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the RPA ground state calculations allow for the inclusion of residual long-range correlations 
necessary to obtain a realistic, ab initio-type description of nuclei. 

On the other side, models based on phenomenological nuclear interactions have been well 
established over the past decades, and have been very accurately tuned to the properties of 
finite nuclei. These models are nowadays successful not only in the region of stable nuclei, 
but also in the description of exotic nuclear structure and collective excitations in nuclei 
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ley of /3-stability, both in the non-relativistic 
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201 ] and relativis- 



221 ]. Particularly interesting is the phenomenon of the pygmy dipole 



resonance in neutron-rich nuclei (PDR), indicating that the loosely bound neutrons might 
coherently oscillate against the approximately isospin-saturated proton-neutron core [23]. 
Very recently it has been shown that for proton rich nuclei in the lower region of the nuclide 
chart, one also could expect the appearance of a low-energy exotic collective mode, i.e. the 



proton PDR, when loosely-bound protons vibrate against the rest of the nucleons [2J|. The 
experimental evidence about low-lying excitations in nuclei towards the drip-lines is still 
rather limited, and actually available only in light nuclei up to the oxygen isotopes [2.5J. 
The present article will summarize the recent progress in studies of low-lying excitations 
towards the proton drip-line within the Relativistic Quasiparticle RPA (RQRPA). 



II. THE HARTREE-FOCK MODEL IN THE UCOM FRAMEWORK 



The UCOM approach aims at an explicit treatment of interaction-induced short-range 
central and tensor correlations in nuclei [llj, |l2|, [l3|, |l4 ] . These correlations are imprinted into 



an uncorrelated many-body state through a state-independent unitary transformation 
defined by the unitary operator C, resulting in a correlated state, \^f) = C \^/). Due to the 
unitarity of the correlation operator, matrix elements of an operator O in correlated many- 
body states are equal to those evaluated using the correlated operator 6 and uncorrelated 
many-body states, i.e. 



(|r| O ]*') = (tf | C f O C \^') = (tf | 6 \^'). 



(1) 



The short-range central and tensor correlations are separately included via the unitary op- 
erators C r and Cn, respectively, and formulated as exponential functions of the two-body 
Hermitian generators g r and g^. The operator form of the generators is motivated by the 
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basic physics of the two types of correlations we are going to describe explicitly. 

The short-range central correlations are induced by the strong short-range repulsion in 
the central part of realistic NN interactions. This repulsive core prevents nucleons in a 
many-body system to approach each other closer than the characteristic size of the core. 
In the two-body density matrix these correlations are revealed through the depletion of 
the probability density for particle distances smaller than the core radius. In order to 
include these correlations into an uncorrelated many-body state, e.g., the Slater determinant 
of the Hartree-Fock approach, we perform a distance-dependent radial shift with respect 
to the relative coordinate of two-particles. The corresponding Hermitian generator g r = 
|[s(r)q r + q r s(r)] contains the radial component of the relative momentum operator, q n 
and a function s(r) which controls the distance- dependence of the shift. Implementation 
of the correlation operator in coordinate representation corresponds to a norm-conserving 
coordinate transformation r i— > R_{r)~ of the relative coordinate. The radial correlation 
function R~(r) and its inverse R+(r) are related to the shift function s(r) 

For the following calculations the Argonne V18 (AV18) potential [q] is used. In practice, 
the correlation functions R+(r) are parameterized and the optimal parameters are deter- 
mined for each spin-isospin channel from an energy minimization in the two-body system 
Q|. Since the (S = 1,T = l)-channel of the AV18 potential is purely repulsive, we employ 
a simple constraint on the range of the correlation function in this channel in order to avoid 
artificial long-range correlation functions. We have checked explicitly that the effect of vari- 
ations of this constraint around the value that ensures the short-range of the correlation 
function, is negligible for all calculations presented here. 

A second, equally important type of correlations is induced by the tensor part of the 
interaction [12]. They entangle the relative spatial orientation of two nucleons with their 
spin orientation. The generator gQ has to describe an angular shift depending on the spin 
orientation. This is achieved by gn = |i?(r)[(<r 1 ■ q n )(cr 2 - r) + (o^ ■ r)(<r 2 - q^)], where 
q n = q— ^q r Q|. As for the central correlations, the tensor correlation functions $(r), which 
control the distance-dependence of the tensor correlator, are parameterized and determined 
by a two-body energy minimization. A characteristic of the tensor part of realistic NN 
interactions and thus of tensor correlations is their long range. We are not aiming at a 
description of long-range tensor correlations by the unitary transformation, since they are 
strongly system-dependent. In contrast to the deuteron, the long-range tensor correlations 
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in heavier nuclei will be largely screened. During the determination of the tensor correlators 
we therefore constrain the range or volume of $(r) given by 

jf T ) = [ drr 2 $(r). (2) 



(S=l T=0) 

In the present study, we use the optimal tensor correlators for I# = 0.07, 0.08, 

and 0.09 fm 3 and investigate the effect on the global properties of collective excitation 
phenomena. No tensor correlator is employed in the (S = 1,T = l)-channel, because 
there the tensor interaction is rather weak. It has been shown within no-core shell model 
calculations that for ij 1 " = 0.09 fm 3 the experimental binding energies for A < 4 are 
reproduced quite well (l4|. 

The correlated operators constructed by the unitary transformation contain irreducible 
contributions, not only of one-body and two-body operators, but also higher n-body parts. 
This cluster expansion is truncated after the two-body level, leading to the so-called two- 
body approximation. In previous studies it has been verified that higher order contributions 
due to central correlations can be neglected in the description of nuclear structure properties 
[l^ . For the tensor correlators the range constraint is important for the validity of the two- 
body approximation. The size of residual three-body and higher order contributions was 
estimated in jl^ . 

Starting from the uncorrelated Hamiltonian for the A-body system, consisting of the 
kinetic energy operator and the bare AV18 potential, the central and tensor correlation 
operators are employed to construct the correlated Hamiltonian in the two-body approxi- 
mation. Therein the one-body contributions come only from the uncorrelated kinetic energy, 
while two-body contributions arise from the correlated kinetic energy and the correlated po- 
tential, which together constitute the low-momentum correlated interaction Vucom HQ. 
The Vucom interaction can directly be employed in the HF model to determine the single- 
particle wave functions and energies. By using the expansion in the harmonic oscillator 
basis, the HF equations are solved in a self-consistent way, with restrictions on the max- 
imal value of the major shell quantum number iV max = 12, and maximal orbital angular 
momentum quantum number Z max = 8. 

In Fig. ^the UCOM-HF single-nucleon spectra are displayed for the case of 40 Ca. The cal- 
culations are based on the correlated Argonne V18 interaction, with the constraint 1^=0.09 
fm 3 for the correlation volume of the tensor correlator, Eq. (j2J). The UCOM-HF energy levels 
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are compared with the HF spectrum obtained with the low-momentum NN potential V^-k. 

with two standard phenomenological interactions in the nonrelativistic (Skyrme) 27| 
and relativistic (NL3) 28] framework, and with experimental levels j^. The spectra ob- 
tained from the HF model based on realistic NN interactions appear distributed too wide 
in energy. In addition, the HF binding energies and radii are too small compared to the 
experimental values |l3| |. 

These deviations can be attributed to several missing pieces in the UCOM-HF description: 
(i) Long-range correlations are not covered by the unitary correlation operators and should 
be described by the model space, i.e. the available many-body states. The independent- 
particle states of the HF approach are clearly not able to do so and one has to go beyond the 
mean-field level. In the next section we are going to include long-range correlations due to 
collective vibrations, by means of RPA. (ii) Three-body forces, either genuine or induced by 
the unitary transformation, generally play a role for the quantitative description of nuclear 
structure. These are not included in the present study. Hence the results presented here 
will provide some information on their importance. 



III. RANDOM-PHASE APPROXIMATION BASED ON THE Vucom 

The HF description of the nuclear ground state is to some extent oversimplified, and 
correlation effects going beyond mean-field should be included. Giant resonances may have 
some influence on the nuclear binding energies ^ | , and it is known that correlations due to 
surface vibrations have a considerable influence on the ground state densities [33]. In this 
section, we will employ an RPA model based on the UCOM Hamiltonian (UCOM-RPA) to 
evaluate the ground state correlations due to collective vibrations and to study the properties 
of such excitations themselves. 

The UCOM-HF single-particle states are used for the construction of the ph configuration 
space for the RPA model. One of the standard approaches to derive the RPA equations is 
the equation of motion method with the quasiboson approximation Q], resulting in the 
eigenvalue problem formulated as a set of coupled equations for the forward and backward 
amplitudes, X^ M and Y^ JM respectively, 

m n \ / x k ' JM \ 

(3) 
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The eigenvalues correspond to RPA excitation energies and the RPA matrices are given 
by, 

= (01 [[^.Hucom] ,A$ + ] 10) (4) 
B J php , h , = [K.Hucom] ,(-1) J - m ^] 10), (5) 

where the operator (A^f) creates (annihilates) a p/t state of angular momentum JM. 

We consistently use the intrinsic Hamiltonian Hucom = T — T cm + Vucom, i-e., the center 
of mass contribution to the kinetic energy is subtracted on the operator level. The Coulomb 
interaction is included explicitly. 

An essential property of the present model is that it is fully self-consistent, i.e. the same 
correlated realistic NN interaction Vucom is used in the HF equations that determine the 
single-particle basis, and in the RPA residual interaction entering the calculation of the 
RPA matrices. This essential property of our model ensures that RPA amplitudes do not 
contain spurious components associated with the center-of-mass translational motion. We 
have verified that the spurious 1~ state is properly decoupled from the physical excitation 
states. We also have examined, for closed-shell nuclei across the nuclide chart, that the 
UCOM-RPA m0 de. essentia!* e, haU sts the isoscalar e„e rgy -wei g Med su m r „!es Q with 
maximal discrepancies of ± 3%. 

In the present study, the correlation energies are evaluated within the UCOM-RPA frame- 
work, 

5i? = -^(2J + l)^^|y; h J | 2 , (6) 

k, J ph 

by using the RPA eigenvalues lu£, and backward-going amplitudes Y jf ' . Both, the natural 
7i = (— 1) J and unnatural parity ix = (— 1) J+1 excitations are included, in the range of 
= ± - 10 ± . The HF binding energies together with RPA correlations due to collective 
excitations are shown in Fig. El for several closed-shell nuclei. The HF binding energies 
and UCOM-RPA correlations are calculated in a consistent way by using the correlated 
AV18 interaction for various ranges of the tensor correlator, constrained by J^ 10 ^=0.07, 0.08, 
and 0.09 fm 3 . In general, the tensor correlator with longer range provides stronger binding 
both on the HF level and when the correlations are taken into account. In comparison 



with the experimental binding energies 



the present model with full implementation of 



RPA correlations seems to favor the tensor correlator with shorter range. However, one 
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should keep in mind that the method used to evaluate the correlation energy, Eq. (0), is 
not free of over-counting and therefore the correlation effects are overestimated. Within 
many-body perturbation theory or configuration interaction calculations, the longer ranged 
correlator provides a very good agreement with experimental binding energies for all nuclei 
in accord with the no-core shell model calculations discussed in 

By examining the electric transition strength, one can also study the properties of the cor- 
related interaction Vucom- This is exemplified in Fig. |3J where the UCOM-RPA strength 
distributions, corresponding to the isoscalar giant monopole resonance (ISGMR) are dis- 
played for the correlated AV18 interaction with different restrictions on the range of the 
tensor correlator, 1^=0.07, 0.08, and 0.09 fm 3 . For the lighter nuclei 16 0, 40 Ca, and 48 Ca, 
the ISGMR is fragmented into two-three peaks, whereas for 90 Zr, 132 Sn, and 208 Pb the IS- 
GMR is strongly collective, resulting essentially in a single peak. For a comparison, the 
monopole response is also calculated in the framework of relativistic RPA based on effective 
Lagrangian with density-dependent meson-nucleon vertex functions, with DD-ME1 inter- 
action (more details are given in Sec. HVj) . In addition, the calculated ISGMR strength 
distributions are compared with the nonrelativistic RPA based on Woods-Saxon potential 



29 



30| and ( 3 He, 3 He) 



and G-matrix formalism |32j, and with experimental data from (a, a) 
scattering One can observe that a decrease of the range of the tensor correlator system- 
atically pushes the transition strength towards lower energies. In particular, by decreasing 
the range of the tensor correlator, i.e. its constraint 1^=0.09 fm 3 towards 0.07 fm 3 , the 
excitation energy of ISGMR lowers by ~4 MeV. This means that by varying the range of 
the tensor correlator, one can effectively control the impact of the missing correlations on 
the transition strength of ISGMR. 



IV. PROTON PYGMY DIPOLE RESONANCE IN THE RELATIVISTIC QRPA 

In this section we discuss recent developments regarding the response of nuclei far from 
/3-stability. One of the major challenges in this region is the understanding of soft modes 
of excitations which involve loosely bound nucleons. In particular, in neutron-rich nuclei, 
nucleons from the neutron skin may giv e rise to a soft low-energy dipole mode known as the 
pygmy dipole resonance (PDR) (3, 0, E3] . The structure of nuclei on the proton-rich side 
is equally important for many aspects of the underlying many-body problem and effective 
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nuclear interactions. In a recent relativistic QRPA study, it has been predicted that in nuclei 
with proton excess one could expect the appearance of a p roton PDR mode where loo sely 
bound protons vibrate against the rest of the nucleons 

HQ- The relativistic QRPA Q 
is formulated in the canonical single-nucleon basis of the relativistic Hartree-Bogoliubov 
(RHB) model and is fully self-consistent. For the interaction in the particle-hole channel 
effective Lagrangians with nonlinear meson self-interactions or density-dependent meson- 
nucleon couplings are used Q|, and pairing correlations are described by the pairing part of 
the finite-range Gogny interaction |2|. The parameters of the effective relativistic interaction 
have been adjusted to properties of nuclear matter and to binding energies, charge radii, and 
differences between neutron and proton radii of spherical nuclei |4QJ]. In the small-amplitude 
limit, the RQRPA equations are derived from the equation of motion for the generalized 



nucleon density 2l|. The RQRPA configuration space is constructed from standard (2qp) 



)airs, but one also needs to include transitions to the unoccupied states from the Dirac sea 

im 

The RQRPA dipole strength distributions for N=20 isotones, displayed in Fig. EJ are 
dominated by the isovector giant dipole resonances (GDR) at ~ 20 MeV excitation energy. 
With the increase of the number of protons, low-lying dipole strength appears in the region 
below the GDR and, for 44 Cr and 46 Fe, a pronounced low-energy peak is found at ~ 10 MeV 
excitation energy. In the lower panel of Fig. 0] we plot the proton and neutron transition 
densities for the peaks at 10.15 MeV in 44 Cr and 9.44 MeV in 46 Fe, and compare them with 
the transition densities of the GDR state at 18.78 MeV in 46 Fe. Obviously the dynamics 
of the two low-energy peaks is very different from that of the isovector GDR: the proton 
and neutron transition densities are in phase in the nuclear interior and there is very small 
contribution from the neutrons in the surface region. By exploring the RQRPA amplitudes, 
we note that, rather than a single proton 2qp excitation, the low-lying states are characterized 
by a superposition of a number of mainly proton 2qp configurations. The low-lying state 
does not belong to statistical El excitations sitting on the tail of the GDR, but represents 
a fundamental mode of excitation: the proton electric pygmy dipole resonance (PDR). 



10 



V. SUMMARY 

In the present study, a fully self-consistent RPA model is constructed in the single- 
nucleon Hartree-Fock basis, by using correlated realistic NN interactions obtained within the 
UCOM framework. It is shown that the Vqcom interaction generates a strongly collective 
ISGMR mode, whose energy is sensitive to the range of the tensor correlator. The UCOM- 
RPA correlations due to collective vibrations provide important contributions to the nuclear 
binding energies. In addition, by employing the fully self-consistent relativistic quasiparticle 
RPA, it is indicated that the nuclei towards the proton drip-line are characterized by the 
appearance of the proton pygmy dipole resonance, i.e. an exotic mode where loosely bound 
protons oscillate against the isospin-saturated proton-neutron core. 
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FIG. 1: The UCOM-HF neutron and proton single particle spectrum for 40 Ca, along with the cor- 
responding spectra from the HF model based on the low-momentum NN potential V\ ow -k ,26], HF 
with SIII Skyrme-type interaction J27 1 , relativistic mean field theory with NL3 effective interaction 
2sl |. and experimental spectrum The UCOM-HF calculations are based on the correlated 

Argonne V18 interaction (/1 10 ^=0.09 fm 3 ). 
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FIG. 2: The Hartree-Fock binding energies for a series of closed-shell nuclei, in comparison with 
the corrected binding energies due to UCOM-RPA correlations from collective vibrations, and 



experimental data |3j]. The binding energies with the full correction (dashed line) and only with 
natural parity excitations (dotted line, n.p.) are separately displayed. Both the Hartree-Fock and 
RPA calculations are based on the correlated Argonne V18 interaction with various ranges of the 
tensor correlator ( 4 10) =0.07, 0.08, and 0.09 fm 3 ). 




FIG. 3: The UCOM-RPA monopole transition strength distributions for the correlated Argonne 
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and 0.09 fm 3 ). The grey lines correspond to the monopole response from the relativistic RPA, 
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ME1 parameterization 
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FIG. 4: The RHB+RQRPA isovector dipole strength distributions in the N=20 isotones, calculated 
with the DD-ME1 effective interaction. For 44 Cr and 46 Fe the proton and neutron transition 
densities for the main peak in the low-energy region are displayed in the lower panel and, for 46 Fe, 
the transition densities for the main GDR peak. 



